Surface spin fluctuations probed with flux noise and coherence 
in Josephson phase qubits 



O 
(N 

> 
O 

(N 



CD 



o 



> 
o 

oo 



X 



Daniel Sank^, R. Barends'^, Radoslaw C. Bialczak^, Yu Chen^, J. Kelly^, M. Lenander^, E. 
Lucero\ Matteo Mariantonil^^ M. Neeley^''', RJ.J. 0'Malley\ A. Vaisencher\ H. Wang^'^^ 
Wenner\ T.C. White\ T. Yamamoto^, Yi Yin\ A.N. Clelandl'^ and John M. Martinisi'^] 
^Department of Physics, University of California, Santa Barbara, California 93106-9530, USA 

^Department of Physics, Zhejiang University, Hangzhou 310027, China 
^ Green Innovation Research Laboratories, NEC Corporation, Tsukuba, Ibaraki 305-8501, Japan 
'^Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, MA 024-20-9108 and 
^California NanoSystems Institute, University of California, Santa Barbara, CA 93106-9530, USA 

(Dated: November 15, 2011) 

We measure the dependence of qubit phase coherence and surface spin induced flux noise on 
inductor loop geometry. While wider inductor traces change neither the flux noise power spectrum 
nor the qubit dephasing time, increased inductance leads to a simultaneous increase in both. Using 
our protocol for measuring low frequency flux noise, we make a direct comparison between the flux 
noise spectrum and qubit phase decay, flnding agreement within 10% of theory. The dependence of 
the measured flux noise on inductor geometry is consistent with a noise source correlation length 
between 6 and 400 /^m. 



Superconducting qubits [l[ are rapidly approaching the 
requirements needed for fault tolerant quantum compu- 
tation, with recent experiments demonstrating a set of 
information storage, logic gates, and coherence improve- 
ments 0-13 • However, in recent demonstrations of both 
multi-qubit gates with phase qubits |6|, |7|] and error cor- 
rection with transmons @, the fidelity of the quantum 
process was limited by individual qubit dephasing times 
T2. This dephasing results from low frequency fluctu- 
ations in the magnetic flux used to tune the devices' 
|0) — |1) transition frequency /iq. Traditionally, flux, 
transmon, and quantronium qubits have addressed this 
problem by idling at "sweet spots" where the first order 
sensitivity of device frequency to bias flux, d/io/d$, is 
zero [Qhllj. This approach has limitations, however, as 
qubits have to be moved from their idle points to partic- 
ipate in multi-qubit quantum gates. Furthermore, this 
approach to the flux noise problem constrains the qubit 
design space, as devices must be carefully engineered to 
have a sweet spot. 

Anomalous low frequency flux noise in superconduct- 
ing devices has a long history in both experiment and 
theory [iMll. Experiments with SQUIDs have found 
that the flux noise power spectral density S'$(/), which 
scales with frequency approximately as 1//, is insensitive 
to device size and materials; in fact no strong depen- 
dence on any device parameter has been found. Recent 
experiments have identified spins on the metal surfaces 
as the noise source 16, but its exact microscopic 
mechanism remains unknown. Experiments with qubits 
have measured the dependence of dephasing on the bias 
point, characterized by dfwjd<^ |l8l - l2d |. and have made 
direct [2l| and indirect [23| measurements of the noise 
spectrum, but no experiment has carefully compared de- 
phasing times with a direct measurement of the low fre- 
quency flux noise, nor has any experiment demonstrated 



a means of improving T2 independently of the bias point. 

In this Letter, we present measurements of dephasing 
times and flux noise in phase qubits designed to improve 
T2. In a circuit with increased inductance, we find an 
improvement in phase coherence at all bias points. We 
introduce a protocol for directly measuring low frequency 
flux noise and use this protocol to make a clear and 
quantitative comparison between qubit dephasing and 
the measured noise spectrum. With this protocol, we 
also report the flrst strong dependence of the flux noise 
on device geometry. 

In order to motivate our redesigned qubits, we first re- 
view dephasing theory. Dephasing is caused by random 
fluctuations in the qubit 's transition frequency /lo, char- 
acterized by a spectral density Sf-^g{f). For a Ramsey 
fringe experiment, the theoretical prediction for the time 
dependent decay of the qubit population probability is 
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where /m is a lower cutoff frequency equal to the in- 
verse of the total experiment time [24j . In phase qubits, 
frequency fluctuations are dominated by flux noise with 
a nearly 1// spectral density, Sg,{f) — S^/ , a«l 
[2l| . resulting in a frequency noise spectral density 
%o(/) = (rf/io/d*)'^;//". Inserting this into Eq.®, 
performing the integral for the case a — 1, and adding 
the contribution from energy loss (Ti), yields 
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The log factor is weakly dependent on t and has 
a numerical value of ^24 in the experimental range 
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Figure 1: (color online) Phase qubit devices used in the exper- 
iments. Panel (a) shows the readout SQUID, inductor coils, 
bias coils, and shunt capacitor; panels (b), (c), and (d) show 
close-ups of the inductors and SQUIDs. (b) Standard qubit 
coils with trace widths w = 1.5 fim and L = 710 pH. (c) 
Wide trace coils with w — 6 /im and L = 720 pH. (d) High 
inductance coils with w = 1.5 /im and L = 1330 pH. 
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Figure 2: (color online) Plot of Ramsey decay time T2 vs. 
nonlinearity A for each type of device. Each point comes from 
a fit to Eq. ((3|. For each bias flux, Ti and A are measured 
in separate experiments. Dashed fit lines are drawn with the 
predicted slope of —3/4 in order to compare magnitudes. The 
fit lines for the standard and wide trace designs fall on top of 
one another. The high inductance devices show an increase 
in T2 by a factor of 1.3. Confidence bounds, typically 2% of 
data values, cannot be seen on this scale. 



10 ns < t < 400 ns |25|. Definmg the dephasmg time T2 
by the equation 
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we find 
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T^cx^r^'M^) =Sr''L[ri^] , (4) 



where L is the inductance of the qubit loop and / is 
the loop current. Evidently T2 can be increased by de- 
creasing the sensitivity dfio/dl, but this incurs a decrease 
in the qubit's nonlinearity A, defined as A = /21 — /lo- 
Nonlinearity is a critical figure of merit, as greater non- 
linearity allows for shorter control pulses and therefore 
more quantum gates during the qubit's lifetime [26'|. Wc 
therefore re-express dfio/dl in terms of A according to 
dfio/dl cx A^/" [Hi, and substitute into Eq. (g]), yielding 



ln(T;) = -Jln(A)+ln(L)-iln(5;) 



K. 



(5) 



Here K depends on other device parameters, such as ca- 
pacitance and junction critical current, which control the 
device's operating frequency range. As we do not wish 
to change the device frequency, K must remain fixed. 
Similarly, as we wish to improve T2 without sacrificing 
nonlinearity, we regard A as a scaling parameter. There 
remain two ways to improve Tj*: raise L or decrease 5*1,. 

Each of these two methods was implemented though 
a redesign of the qubit inductor loop. The devices are 
pictured in Fig.[TJ The standard phase qubit (27l . [28| is 
shown in Fig.[T^, with a close-up of the inductor coils in 



Fig.Hb. Coils of the first redesign are shown in Fig.[T}:. 
Theories attributing flux noise to localized, uncorrelated 
magnetic moments of unpaired electron spins on the 
metal surfaces predict that scales a.s R/W where R is 
the radius of the metal loop and W is its width 1^ 21 1 . 
Therefore, to lower S^, we increased the trace widths 
from the standard 1.5 /im to 6.0 /im. The second re- 
design is shown in Fig.[TJi. In this device, we increase the 
inductance of the loop by adding turns, increasing from 
710 pH to 1330 pH. Increased inductance should reduce 
noise currents driven by noise flux and lead to better Tj* 
as per Eq. ([5]). 

In the first set of experiments we measured Tj* sev- 
eral flux bias points that varied A. The exponential en- 
ergy decay time Ti and nonlinearity A are first measured 
in separate experiments (29l . [s^- We then measure a 
Ramsey decay curve and fit the data to Eq. ^ so that 
the contribution from Ti is appropriately removed in ob- 
taining the Gaussian decay constant Tj*. 

The results shown in Fig.[5]reveal that the dependence 
of the extracted values of Tj* on A is consistent with 
the expected -3/4 power law. This indicates that Tj* 
scales inversely with the qubit frequency sensitivity, as 
predicted by Eq. ^ . The wide trace devices show no 
difference in from the standard design over the mea- 
sured range of operating points. Because the inductance 
of the wide trace design is equal to the inductance of the 
standard design, the lack of change in Tj* suggests that 
was not lowered by the increase in trace width, con- 
sistent with previous results in SQUIDs |12j]. The high 
inductance device showed an increase in by 30% rel- 
ative to the standard design at all bias points. This is 
our first main result: increased inductance improves the 
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Figure 3: (color online) The Ramsey tomography oscilloscope 
(RTO) protocol, (a) Four Ramsey sequences are performed, 
two with the final n/2 pulse about the ±X axes, and two with 
the final n/2 pulse about the ±Y axes. The precession time 
for both sequences is fixed at r = 100ns. (b) The combined 
X axis and Y axis sequences constitute tomography of the 
state in the x-y plane after precession by r, thereby measur- 
ing the angle ip traversed by the state, (c) Several hundred 
subsequent X and Y measurements are grouped together to 
get an average value for if. The entire process is repeated 
once per second for several hours to build up a time series 
Sf{t) = Sip{t)/27rT. 



qubit's dephasiiig time independently of the bias point, 
and without sacrificing any other figure of merit. 

The disparity between the l.SOx measured improve- 
ment and the 1.87x improvement predicted by Eq. ([SJ, 
based on the designed increase in L, suggests that the 
high inductance device had an increase in flux noise. Ac- 
cording to Eq. ^ the relative improvement in T2 is ex- 
pected to be 



(6) 



For the measured Tj* increase of 1.30 and the designed 
inductance increase of 1.87, Eq. ([6|) predicts an increase 

in noise amplitude of ^ /S^ — 1.4. 

We next test these predictions for the relative noise 
changes against a direct measurement of the flux noise. 
The Ramsey fringe measurement is sensitive to the time 
dependent fiuctuations of the qubit frequency averaged 
over the entire data acquisition time, and therefore mea- 
sures integrals of the noise spectral density. To directly 
probe <S'$(/), one can instead track the qubit frequency 
in real time and then produce a spectrum with Fourier 
transform methods. A previous experiment with phase 
qubits implemented this idea by measuring the time de- 
pendent fluctuations of the position of the qubit reso- 
nance peaks |2l| . The present experiment refines this 
idea: instead of measuring /lo spectroscopically, we use 
the free precession of the qubit state in the equator of the 
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Figure 4: (color online) Power spectra for each inductor de- 
sign, (a) Typical spectra for each type of device, (b) Spectral 
power integrated over the range 0.001 Hz to 0.1 Hz. In (a) we 
plot only one spectrum for each device design for clarity in 
the plot. In (b) we show several curves, measured at several 
bias points, in devices all fabricated on the same chip. 



Bloch sphere. The angle traversed by the state in a fixed 
time T is ip = 27r/ioT. By clioosing r and measuring ip 
through partial tomography we obtain /iq. Repetition of 
this measurement for several hours produces a time se- 
ries fio(t). This measurement protocol, which we call the 
Ramsey tomography oscilloscope (RTO), is illustrated in 
Fig.m The RTO requires Httle caHbration and can be 
implemented in any qubit system. 

The time series fio{t) is converted to a periodogram 
using the discrete Fourier transform, and the result is 
appropriately normalized to produce the power spectral 
density of frequency noise Sf-^g{f) The sensitivity 

dfia/d^ is measured separately, and the flux noise spec- 
tral density is obtained from the frequency noise spectral 
density according to S'$(/) — {dfio/d^)^^Sf^g{f)- 
suits are shown in Fig. S) Note in Fig.|4^ that the spec- 
tra contain little statistical noise, in contrast to the data 
of Ref. [21] where the measured power spectra wandered 
by almost an order of magnitude within each decade of 
frequency. 

The power spectra for each type of device scale as f~^'^ 
over the measured band and have extrapolated ampli- 
tudes at 1 Hz between 3.5 and 5.0 /i^o/VSz. The spectra 
in Fig. 2^ directly show that the high inductance device 
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has more flux noise than the other two devices. Still, 
a clearer picture is obtained if we integrate the power 
spectral density to produce a plot of cumulative power. 
To avoid complications from aliasing and spectral leak- 
age from DC, we omit the upper and lower edges of the 
measured frequency band, integrating only over the range 
10"'^ to 10""'^ Hz where the measurement is the least sus- 
ceptible to these effects |31|. From the integrated power 
shown in Fig.|4jD, it is clear that the wide trace and nor- 
mal designs have the same level of noise, while the high 
inductance device shows a 1.7x increase in noise power. 

For the wide trace devices, the lack of change in noise 
level relative to the standard devices is consistent with 
the results of the Ramsey fringe experiments in which 
these two designs had the same values of T2 ■ For the high 
inductance device, the \/1.7 = 1.3-fold increase in noise 
amplitude found in the RTO is in approximate agreement 
with the results from the Ramsey fringe experiments and 
Eq. 16]), which predict an increase of 1.4. 

We have also used the RTO to measure cross correla- 
tion in the noises of two phase qubits separated by 500 
/im on the same chip. We find that the measured correla- 
tion is no greater than that found for two independently 
simulated 1// noise signals j25| . 

The scaling of the noise power with inductance, and 
the lack of scaling with trace aspect ratio, have impli- 
cations for models of the flux noise. The lack of change 
with aspect ratio is incompatible with models predict- 
ing oc R/W, such as the independent surface spins 
proposed in Ref . [2l[ , or single electrons interacting via 
the superconducting condensate as proposed in Ref. 
This may indicate a correlation length larger than the 
trace width. If the noise sources were correlated over the 
entire length of the coil, we would expect the noise power 
to scale quadratically with the number of coil turns, while 
for uncorrelated sources it would scale linearly. The high 
inductance devices had 50% greater number of turns than 
the standard devices, and we found a noise power increase 
of 70%, closer to the uncorrelated case. The near linear 
scaling is also evidence that the noise sources produce 
local fields; non-local fields would allow each source to 
couple spin into multiple turns of the coil, making the 
scaling of noise power with number of turns super-linear. 
A possible interpretation of all of these results is that the 
noise sources are localized, but have a correlation length 
that is greater than the 6 /im trace width and smaller 
than the 400 fj,m. total length of the high inductance coil. 

Thus far we have discussed the relative dephasing times 
and flux noise levels. It remains to check whether the 
the absolute measured noise levels accounts for the ob- 
served dephasing times. This is done by comparing the 
value of g*^^'^^'^y extracted from a flt to the Ramsey data 
with the value found directly using the RTO. Note 

that previous experiments with qubits have been unable 
to make this comparison because an accurate measure- 
ment of the flux noise was unavailable. As described by 



Eq. ID) and the surrounding discussion, the exact form 
of the Ramsey fit function depends on the scaling power 
a. In particular, it would be inappropriate to fit the 
Ramsey data to Eq. ([2]) ; that equation was derived under 
the assumption that a is equal to 1, whereas we found 
in the RTO that a is closer to 1.1. To take the scaling 
power into account properly, we numerically evaluate the 
integral in Eq. JT]) for a — 1.1, and use the resulting fit 
function to extract S^'^"^^'^^ . For all devices and bias 
points we find that the values of g^^™^°y extracted in 
this way agree with 5'!^'^° to within 10%. We empha- 
size that properly accounting for the measured slope of 
the noise spectral density when fitting the Ramsey curves 
is essential in obtaining agreement between the Ramsey 

and RTO experiments. For example, if we take a = 1 we 
find 5;^"'"'''='' « 4 5™^ 

In conclusion, we have found that increased loop induc- 
tance is a viable means to improve superconducting qubit 
phase coherence independently of bias point. We have 
introduced the RTO protocol for directly measuring fiux 
noise, and used this protocol to make a direct comparison 
between measured fiux noise and qubit dephasing times. 
There we find that the noise level extracted from Ramsey 
decay agrees with the direct measurement only when the 
exact slope of the noise spectrum is considered. From an- 
other point of view, this sensitivity means that Ramsey 
decay times can only be accurately predicted when the 
scaling power of the noise is known. Using the RTO we 
compared the fiux noise in devices with widened traces 
and increased number of inductor turns. With no change 
found in the wide traces, but a nearly linear increase in 
noise with number of turns, we conclude that the noise 
sources may be correlated over distances greater than the 
6 /xm trace width, but smaller than the 400 /Ltm length 
of the inductor coils. Clearly, the correlation length of 
the fiux noise sources is a key parameter that should be 
studied in further experiments. 
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Here we present experimental details on the Ramsey Tomography Oscilloscope (RTO) protocol 
and details of the calculations used to extract the flux noise magnitude from Ramsey decay data. 



RAMSEY TOMOGRAPHY OSCILLOSCOPE 

Here we describe the data processing procedure for the 
Ramsey Tomography Oscilloscope (RTO). We found that 
careful signal processing was important in reducing sta- 
tistical noise in the power spectra generated by the RTO. 
The bandwidth of the RTO measurement is set funda- 
mentally by the rate at which the qubit can be measured 
and reset. In our case this would allow ideally 10,000 
quantum measurements per second. With our current 
asynchronous control software this limit could not be 
reached while simultaneously tracking the time at which 
each measurement occurred. Maximum data rate with 
accurate time stamping was achieved with 2,400 quan- 
tum measurements per second, 600 of each of the four 
tomography sequences. Averaging the 600 measurements 
together produced one frequency measurement per sec- 
ond. This set the bandwidth of the experiment to be 0.5 
Hz due to the Nyquist criterion. 

Data was typically acquired for eight to ten hours, 
yielding between 28,000 and 36,000 points in the time 
series. Power spectra are computed as follows. First, the 
time series is divided into four or five non-overlapping 
sections. We compute the power spectrum of each sec- 
tion separately and average them together at the end of 
the procedure. To eliminate uncorrelated quantum mea- 
surement shot noise, we use an interleaving procedure on 
each section. Each section is split into two interleaved 
time series, fi(n) and f2in) {n is the discrete time in- 
dex). These series are multiplied by Hann windows and 
the discrete Fourier transforms Fi (k) and F2 {k) are com- 
puted (fc is a frequency bin index). We form the product 
Fi{k)F2 {k), average neighboring bins together using a 
Gaussian weight function with full width at half maxi- 
mum (FWHM) of 20 bins, and take the magnitude to 
obtain the periodogram P{k). Next, the periodogram 
is multiplied by a factor 1/0.375 to correct for the loss 
of incoherent (noise) power caused by application of the 
Hann window 0] . The periodogram is then smoothed by 
averaging neighboring frequency bins with a Gaussian 
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Figure 1: Cross spectra, (a) Cross spectrum measured using 
the RTO. (b) Cross spectrum computed from two indepen- 
dently simulated 1/f noise signals. 



weight function with a variable FWHM scaling quadrat- 
ically from 1 bin at the low end of the frequency band 
to 20 bins at the high end. The power spectrum S{f) is 
then computed from the periodogram according to 

2T 

where T is the total length of time represented by the sec- 
tion of the time series, and N is the number of points in 



2 



45r 
40- 
35- 

£25- 

£20- 
c 

-15- 
10- 



a= 1.0 
a= 1.07 
a= 1.15 



■ ■ 



» I ; ♦ i 1 1 1 1 



10 

X 



15 



20 



is 



400' 
350- 
300- 
250- 
*~h200- 
150- 
100- 
50- 



A A A A A A A A 4 A A A A A 



a= 1.0 
a= 1.07 
a= 1.15 



A AA i 



50 100 150 200 250 300 350 400 

t [ns] 



Figure 2: (color online) The integrand of /. The curves are 
well behaved over the entire integration range. 



the section. Finally, spectra generated from each section 
are averaged together. 



CROSS CORRELATION 

In order to check that the flux noise we measured was 
generated locally to each device, we used the RTO to 
measure cross correlation of the noise signals generated in 
two devices separated by 500 /im on the same chip. Time 
series of the two devices' resonance frequencies were mea- 
sured using the RTO, and the cross correlation was com- 
puted. Results are shown in Fig. [1] Although there are 
frequencies at which the cross correlation amplitude is as 
high as 0.3, this must be compared against the cross cor- 
relation computed for two independently simulated noise 
signals. We find that the cross correlation of two inde- 
pendently simulated 1// noise signals show very similar 
peak structure to the data, indicating that the noises 
within the two qubits are no more correlated than in- 
dependent noise. This result agrees with the finding in 
Ref. [3| , where it was inferred from quantum state tomog- 
raphy performed on two coupled qubits that dephasing in 
each qubit was uncorrelated. We note that the absence 
of a low frequency roll-off in the RTO data indicates that 
the low frequency flux noise is correlated on time scales 
exceeding the length of data acquisition. For this rea- 
son it is unsurprising that residual cross-correlation was 
found in both the data and the simulation. 



COMPARISON OF RTO AND RAMSEY DECAY 

We wish to fit our Ramsey decay data to the theoretical 
curve given by Eq. (1) in the main text, which, for the 



Figure 3: (color online) The integral / evaluated versus t for 
several values of a. Note the strong sensitivity to q; as a 
goes from 1.0 to 1.15, a 15% change, the integral increases by 
a factor of ~10. 



case where the flux noise is S'$(/) ~ S^/ f", is 



p{t) = exp 



(27r)2 fdf^ 



sin(7rz)^ dz 
(7rz)2 z" 



(2) 

Here sinc(x) = sin(^)/a;, /m ~ 1 hour and t is in the 
range to 400 ns [4|. In order to do this we need to 
evaluate the integral 



sin(7rz)^ dz 
(7rz)2 



(3) 



We compute the integral numerically. Since f^^t is on the 
order of 10""'^^, the lower limit of integration is a very 
small positive number and the integrand is diverging at 
the lower limit. On the other hand, the integrand oscil- 
lates for z > 1. The integral is therefore unfit for numer- 
ical analysis in its current form as it has both divergent 
and oscillatory behavior. The problem is mitigated by 
the change of variables x = — ln(z) which yields 



-ln(/„t) 



sin(7re 



dx 



(4) 



The integrand is now well conditioned over the whole in- 
tegration range. Plots of this integrand for several values 
of a are shown in Fig. (5] Note that an upper cutoff in 
the frequency integral would translate to a lower cutoff 
in the integral over x. Because of the logarithmic scale 
combined with the very small value of the integrand for 
values of a; < 5, ignoring a possible upper cutoff greater 
than 1 MHz incurs negligible error. 

We perform the integral / for 50 values of t in the ex- 
perimental range to 400 ns, and for several values of 
a near 1. Results of the integration as a function of t 
are shown in Fig.[3l We also show / as a function of a 
for two fixed values of t in Fig.|31 From these curves we 
construct interpolating functions and use them to fit our 
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Figure 4: (color online) The integral / evaluated as a function 
of a for several values of t. Note the strong dependence on a. 



measured Ramsey decay data to Eq. Note particu- 
larly in Fig. [4] the strong dependence of the noise integral 
on a. It is because of this strong dependence that we are 
able to accurately determine which value of a gives the 



best agreement with the power spectra measured directly 
using the RTO, as described in the main text. 
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The lower cutoff frequency /m = 1/hour arises from the 
manner in which we acquire the data. Because we use a 
projective quantum measurement to read the state of the 
qubit, we must repeat the experiment for each value of t 
many times to reduce statistical noise. Rather than aver- 
age each point in t sequentially, we spread the averaging 
of each point over the entire trace acquisition period. This 
results in better averaging of the noise signal as explained 
in Ref. Q. 



